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Utilization of Satellite Data for 
Inventorying Prairie Ponds 
and Lakes 


LANDSAT>1 data were used to discriminate ponds and lakes for 
waterfowl management. 


Introduction 

T he primary breeding areas of North 
American waterfowl are the Dakotas, 
the south* portions of the prairie prov- 
inces, nort. v. estem Canada, and parts of 
Alaska. These areas are annually the sites of 
systematic surveys conducted by the U.S. 


lished sampling transects. Numbers and dis- 
tributions of natural ponds and lakes as tabu- 
lated during these flights are an important 
criteria used in assessing annual waterfowl 
production. In order to explore prcxiedures 
which could enhance FWS capabilities for 
monitoring annual and seasonal changes in 


Abstract: By using data acquired by LANDSAT-1 (formerly ERTS- 
1), studies were conducted in extracting information necessary for 
formulating management decisions relating to migratory waterfowl. 
Management decisions are based in part on an assessment of habitat 
characteristics, specifically numbers, distribution, and quality of 
ponds and lakes in the prime breeding range. This paper reports on a 
study concerned with mapping open surface water features in the 
glaciated prairies. Emphasis was placed on the recognition of these 
features based upon water s uniquely low radiance in a single near- 
infrared waveband. The results of this recognition were thematic 
maps and statistics relating to open surface water. In a related ef- 
fort, the added information content of multiple spectral wavebands 
was used for discriminating surface water at a level of detail finer 
than the virtual resolution of the data. The basic theory of this 
technique and some preliminary results are described. 


Fish and Wildlife Service (FWS) in coopera- 
tion with the Canadian Wildlife Service and 
various states and provinces. The surveys 
serve to aid management decisions relating 
to annual hunting regulations and research 
needs by providing an estimate of popula- 
tion size and annual reproduction. These 
surveys rely chiefly upon observations made 
from low-flying light aircraft tra\eliiig estab- 


water conditions, an evaluation of 
LANDSAT-1 (formerly ERTS-1) sensors was 
conducted (Work 1974; Work etal., in press). 
Objectives of that investigation were to map 
and tabulate statistics on surface water con- 
ditions and to determine changes ip wetness 
between spring and summer during a two- 
year period in a glaciated prairie region lo- 
cated in east-central North Dakota. This 
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paper reports on the techniques used to de- 
tect and characterize open surface water fea- 
tures. 

PuLiMiNARY Considerations 

Level thresholding of a radiation signal in 
a near-infrared waveband is a reliable and 
simple technique for delineating surface wa- 
ter. This technique is effective because at 
near-infrared wavelengths the apparent 
radiance of water is usually uniform and 
lower than for other terrain objects. Thus, 
using an appropriate near-infrared 
waveband, water may be delineated by ac- 
cepting scene points with low radiance val- 
ues (classified as water) while rejecting all 
values above a certain threshold (non-w’ater). 
In order to appreciate the use of and limita- 
tions to this technique a brief discussion fol- 
lows. 

The apparent radiance of a body of water 
is the result of (1) reflections at the air-water 
interface, (2) reflections from particulate 
matter suspended in the water volume, and 
(3) reflections from the bottom. Because the 
field-of-view of the L.\NDS.\T scanner is 
limited to near vertical observations (nomi- 
nally to within 5.78 degrees of satellite nadir) 
and because w'ater surQices reflect specular- 
ly, radiation reflected by water to the scan- 
ner must emanate from a sky position near 
the zenith. Given the northerly latitudes 
which characterize the glaciated prairies, the 
LANDS AT scanner does not view water- 
reflected direct solar radiation (i.e., the 
ground specular point is at a considerable 
distance outside the ficld-of-view of the 
scanner).* This leaves only that fraction of 

* Strong (1973), through the combined use of 
imagery collected over ocean surfaces from 
LANDSAT-1 and NOAA-2 satellites, found that 
when the sun elevation exceeds 55 degrees, the 
LANDS AT imagery was subjected to considerable 
contamination by sunlight even though the specu- 
lar point was nearly 550 km from nadir. Stnimpf 
and Strong (1974) termed this condition “diffuse 
glitter." The condition is due to the sea state or the 
slope of wind driven surface waves. In North 
Dakota the maximum solar elevation angle during 
the time of the LANDSAT-l overflight (approxi- 
mately 1015 local sun time) was 59.5 degrees 
along the southern boundary of the state at the 
summer solstice. This would indicate that water 
bodies at these latitudes w’ere capable of produc- 
ing diffuse glitter for a limited period of time cen- 
tered on 22 June. However, surface conditions on 
even the largest water bodies found in the North 
Dakota prairies were considerably smoother than 
those found in the oceans. No evidence of diffuse 
glitter was observed during the course of this 
study. 


diffuse skylight which emanates from a 
near-zenith sky location to impinge ui>on the 
water surface and thence to be reflected to 
the scanner. In relative magnitude, however, 
diffuse skylight is much weaker than direct 
solar radiation, especially in the near- 
infrared and under clear sky conditions 
when optimal satellite observations are pos- 
sible. McDowell (1974) illustrated the mag- 
nitude and spectral differences between dif- 
fuse skylight and direct solar radiation (Fig- 
ure 1). In addition, the reflected skylight 
component is further diminished because 
water surfac'es are a uniformly weak reflector 
of radiation which impinges at any but very 
oblique angles (Figure 2). 

In considering reflections emanating from 
particulates within the water volume and 
from the bottom surface, water’s absorptivity 
must be considered. In the near-infrared, 
that fraction of radiation which penetrates 
the air-water interface is largely absorbed, 
the extent of absorption being dependent 
upon the wavelength and the length of the 
water path. This situation is shown quantita- 
tively in Figure 3, which illustrates the 
spectral transmission of pure water for a va- 
riety of path lengths. Consequently, a sensor 
viewing a water body in a near-infrared band 
receives little or no radiation that may have 
been reflected by the bottom or volume sus- 
pended particulates. 

The relative utility of different near- 
infrared bands for detecting water should 
not be predicated purely on the basis of 
longer wavelengths. For example, the use of 
a waveband in the 2.0- to 2.6-/xm atmos- 
pheric window is not optimal because of the 
decreasing amount of solar radiation at these 
wavelengths. It must be remembered that 
most terrestrial objects are relatively strong 
diffuse reflectors and that with adequate 
solar illumination such targets will contrast 
sharply with darker surface water features. 
An ideal waveband for delineating surface 
water is the 1.5- to 1.8-^m spectral interx’al 
(Work and Thomson 1974; VVork 1974). Suf- 
ficient solar illumination is available in this 
waveband to illuminate background objects. 
By u.sing this waveband, it is also possible in 
many instances to distinguish water which is 
occluded by aquatic plants. These 
capabilities arc of particular importance for 
lower-altitude aircraft data in which case it is 
both possible and desirable to identify small 
ponds, sheet water, and shoreline patterns of 
larger ponds and lakes. Analyzing 
LANDSAT-1 scanner data, however, we 
found band-7 (0.8 to 1.1 fim) to be satisfac- 
tory except for an occasional omission error 


SATELLITE DATA FOR INVENTORYING PONDS & LAKES 


687 



Fic. 1. Scene irradiance components for a clear day. (After 
McDowell. 1974.) Spectral irradiance levels attributable to 
skylight alone may be significantly different and variable widi 
wavelength due to variations in atmospheric haze. The condi- 
tions shown were recorded on an exceptionally clear day in 
New Mexico on 30 October 1970 at a sun elevation angle of 43 
degrees. 


in the case of shallow ponds carrying high 
sediment loads. 

Methods 

Our delineation of surface water was ac- 
complished by using a high-speed digital 
computer which examined each scene pixel 
(picture element) and classified the pixel as 
either water or nnn-water. In training the 
computer for this automated classification. 


data samples were examined in order to de- 
termine typical water radiance levels and 
the radiance levels of other scene materials 
also known to have low radiance characteris- 
tics. Experience has shown that various 
scene objects exhibit radiances that may ap- 
proach the low radiance of water depending 
upon the specific near-infrared waveband 
under consideration, the geographic locale, 
and the phenologic state of scene objects. In 
eastern North Dakota, dark prairie soils. 



ANOLI OP INCmtNCI MIA8UR10 
ritOM NORMAL 

Fig. 2. Variation of reflecting power of the 
air-water interface as a function of the angle 
of incidence (Natural or unpolarized light of 
wavelength 0.589 iim). 



Fig. 3. Spectral transmittance of pure water 
for different path lengths. (Plotted after data 
from Sverdrup et al., 1942.) 
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MoIIisols (formely referred to as Cher* 
nozems), have consistently approached the 
low radiance values of water in tlie 0.8- to 
waveband. In such a situation, sam- 
pled water and soil radiance values ."ere 
examined in order to discern relative differ- 
ences. Typical histograms or frequency dis- 
tributions are shown in Figure 4. In these 
data, both water and soil had distributions 
which were displaced from each other al- 
though slight overlap occurred in the tail re- 
gions. We have found that shallow water fea- 
tures are often represented within the lead- 
ing tail of such a water distribution; as a re- 
sult, die threshold or decision boundary was 
located nearer to the soil’s density peak in 
order to group the shallow water situations 
with other surface water categories. 

Results 

The intent of the recognition process was 
to monitor changes in surface water condi- 
tions for a seasonal (May-to-July) and an an- 
nual period (July-to-July). Data gathered by 
LANDSAT-1 on 31 July 1972, 14 May 1973, 
and 7 July 1973 were analyzed. The same 
area comprising 3311 km* (1278 mi*) was ob- 
served on each of these dates. The study area 
included portions of two different physio- 
graphic regions: a coteau or moraine feature 
created by stagnation ice, and a drift plain or 
low relief feature of numerous ground 
moraines. The drift plain inherently has 


fewer potholes or basin features and, be- 
cause of its low relief, has been subjected to 
numerous wetland drainage projects. This 
difference in wetland occurrence warranted 
a stratification of the numerical results based 
on these physiographic variations. 

A computer-generated tliematic map iden- 
tifying surface water in a portion of the study 
area located almost wholly within the coteau 
physiographic region was produced (Figure 
5). Although a digital map of the type illus- 
trated graphically portrays water conditions, 
it does not represent an efficient form of data 
management. Thematic maps as such, and 



Fic. 4. Histograms used to locate a deci- 
sion or threshold boundarv for separating 
water from non-water terrain features. 



Fig. 5. Digital water recognition map from ERTS observation 1008-16594 of 
31 July 1972. The map, obtained by thresholding channel MSS-7, depicts an 
area northwest of Jamestown, North Dakota. 
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the required manual interpretation, are not 
feasible for wide area surv'eys on a routine 
basis. For purposes of analysis, the data gen- 
erally has greater value if it is presented in 
terms of statistical tables. Figure 6 illustrates 
a computer-generated statistical tabulation 
of ponds and lakes observed on 7 July 1973 
in that portion of the study area lying in the 
coteau physiographic region. The upper list- 
ing itemizes all ponds recognized in the 
scene and provides the location of each pond 
based on a coordinate system derived from 


the sensor’s scanning geometry and also on a 
geographic coordinate system. A size fre- 
quency distribution of ponds in the entire 
scene is presented in the lower tabulation. 

A graphical summary of pond frequency 
for two dates illustrates the pond size dis- 
tribution for the coteau portion of the study 
area (Figure 7). The two curves show that a 
decrease in pond numbers occurred in the 
1973 May-to-July interval. Information on 
changes in pond numbers over this seasonal 
period is an important data input to models 



ERTS-1 at 1654 GMT on 7 July 1973. 
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current'y used for predicting waterfowl pro- 
duction (Geis et ai, 1969). 

In examining these data it should be rec- 
ognized that the pond sizes listed by the 
computer should in practice be termed “ap- 
parent size.” This is because each pixel of 
data was examined and determined to be 
either totally water or not water. Many pixels 
lying on the perimeters of ponds and lakes 
undoubtedly contained some unrecognized 
and untabulated water. This caused the sur- 
face areas of virtually all water features to be 
underestimated. In terms of percentage, the 
error were greater for the smaller ponds and 
for those of irregular shape (i.e., those having 
a high ratio of perimeter length to area). The 
very small ponds, of course, would not be 
recognized. Generally a pond must have 
been at least 0.4 hectare (I.O acre) in size to 
be recognized. Recognition of a 0.4 hectare 
pond was dependent upon whether the pond 
was wholly included in one digital sample 
(i.e., within a pixel) or fractionally distrib- 
uted over several pixels. In general, it is 
problematic whether ponds in the 0.4- to 
I.6-hectare (I.O- to 4-acre) size class were 
recognized. Above 1.6 hectares, ponds were 
nearly always recognized but not necessarily 
recognized at their full areal extent. 

Improvements in Spatial Resolution 

Because prairie ponds are frequently 
smaller than 0.4 hectare (1.0 acre), it was ap- 
parent that many water bodies were not 


tabulated. We , therefore, undertook a radi- 
cally different recognition technique which 
offered promise for improving the data's res- 
olution limit. The remainder of this paper 
describes that approach and discusses some 
preliminary lesults. 

When the instantaneous-field-of-view 
(IFOV) of a scanner is large with respect to 
the scene objects being scanned, a single res- 
olution cell may contain a mixture of mate- 
rials. In attempting to impntve upon the da- 
ta’s resolution capabilities, we have utilized 
the added information content of multiple 
spectral wavebands to estimate the propor- 
tions or fractions of materials which fill the 
IFOV. The technique termed “proportion 
estimation" or “mixtures extimatioii” was 
first outlined by Horwitz et al. (1971) and 
further described by Nalepka et al. (1972). 
Before this study, the application of this 
technique was largely developmental in na- 
ture, and its use in this and another 
LANDSAT-1 study (Malila and Nalepka, 
1973 and 1974) must be considered to be 
among the first attempts to test its applicabil- 
ity in an operational context. 

A discussion of proportion estimation 
theory is beyond ihe scope of this paper, 
however the essenci* of the technique can be 
described in geomet'ic terms. Assume that a 
data set comprised of two spectral channels, 
X| and X}, contains three pure and unique 
materials — A, B, and C. This situation is de- 
picted in Figure 8 where the signature 



Fic. 7. Summary of size distribution of ponds in the coteau stratum 
for a seasonal period. Where the pond size increments are greater 
than one-acre, the data have been averaged over the increment. 


SATELLITE DATA FOR INVENTORYING PONDS & LAKES 


691 


means for the three materials arc shown in 
two-dimensional signal space. In the non- 
degenerate case, each pure signature 
is a distinct vertex of the closed geometric 
figure or signature triangle. If an unknown 
scene element (IFOV) consists of a mixture 
of all three materials, the signature gener- 
ated by this unknown element, X, lies within 
the simplex. An estimate of the proportion of 
each pure material constituting the unknown 
element is obtained by drawing a line from a 
vertex through the signal to be classified to 
the opposite leg of the simplex. The fraction 
of the line between the crossover point and 
opposite leg defines the proportion of the 
corresponding vertex material in the un- 
known. For the case illustrated in Figure 8. 
the unknown happens to lie at the centroid 
of the triangle, and its composition would be in 
the ratio of 1/3, 1/3, and 1/3 of materials A. B. 
and C respectively. A case requiring special 
geometric interpretation is shown in Figure 
9. In this instance the unknown, Z, lies out- 
side or on the edge of the signature simplex. 
The unknown is determined to be com- 
prised of only materials A and C in the in- 



Fic. 8. Geometric interpretation of means 
of signature mixtures. In the case illus- 
trated, the unknown, X, is a mixture of three 
pure materials - A, B, and C - which form 
the vertices of the sigruture simplex. 



Fic. 9. Geometric interpretation of esti- 
mate for a special case. The unknown. Z, 
lying outside the signature simplex is a 
mixture of materials A and C. 


verse ratio by which the simplex triangle's 
leg A-C is divided by a line drawn from Z 
orthogonally to that leg. If the unknown is 
quite distant from the signature simplex (de- 
scril>cd in terms of a distance), the al- 
gorithm designates the unknown as an alien 
object or an object composed of none of the 
simplex materials. 

Although tlie above description has been 
limited to tliree pure and unique materials in 
two-dimensional space, the concept is easily 
expanded to situations where many object 
materials exist in a spectral hyperspace. In 
applying the algorithm, however, it is neces- 
sary that two operational constraints he ob- 
served: (1) at least n-1 spectral channels of 
information are required to satisfactorily es- 
timate mixtures of n-object maL'rials, and (2) 
the signatures for the materials in a mixture 
must not be similar and no one signature 
must come close to the weighted average of 
the other signatures. If either of these latter 
conditions occur, the signature simplex is 
said to be degenerate, and the estimate of 
proportions may be poor or invalid. 

In our study, the intent was to delineate 
water in a mixture of scene materials, 
thereby making it possible to improve the 
size estimates of larger ponds and, more sig- 
nificantly, to recognize small ponds which 
otherwise would not have been detected. 
Examining the data collected in the July 
1973 time frame, it appeared that only two of 
the four LANDSAT wavebands were unique 
in representing the majorit>’ of scene mate- 
rials present in eastern Nortli Dakota. Bands 
5 (0.6 to 0.7 fim) and 7 (0.8 to 1.1 fim) ap- 
peared to be best for distinguishing most 
scene materials while hands 4 (0.5 to 0.6 fxm) 
and 6 (0.7 to 0.8 fim) were respectively re- 
dundant. This being the case, it appeared 
unlikely that the proportion estimation al- 
gorithm could function in using more than 
three signatures. As a result, we chose for 
processing signatures representing water, 
bare soils, and green vegetation. Using these 
signatures, the computational algorithm was 
applied to an area of 287 Ion* ( 1 10 mi*). The 
resultant computer output was a set of water 
proportions for each scene pixel. A water 
recognition map generated from this output 
is shown in Figure 10 and, for comparison, a 
map generated by using the single 
waveband thresholding algorithm is shown 
in Figure 11. In the proportion estimation 
map, the symbol density is related to the 
proportion of water estimated for that pixel. 
In order for the map accurately to portray the 
scene, certain percentage or acceptance 
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limits were determined for the output of the 
algorithm. For example, It seemed ai)i>ro- 
priate to cx)imt pixel values of per cent and 
above as totally water. This procedure 
tended to account for the lil^elihood that a 
value close to the signature mean (i.e., close 
in terms of its probability of ireing that mate- 
rial) nray in fact have been a pure sample 
related to that mean. Similarly, pixels show- 
ing less than 30 per cent water were as- 
sumed to contain no water at all. 

In general, a detailed comparison of the 
classillcation maps and related imagcr>' indi- 
cated that proportion estimation signifi- 
cantly improved ixmd shu)X' definition. Fish 
Lake and several nearby lakes slurnm in F'ig- 
ures 10 and 11 illustrate thi^ fact. .Alkali 
Lake (F'igurc 10), although recognized, was 
not accurately delineated. As implied by its 
name, the lake apparently was high in dis- 
solved and precipitated solids and as such 
was a lacuutrinc feature fur which the water 
signature utilized was not representative. An 
inspection of this same lake in aerial photog- 
raphy indicated a color hue similar to the 
highly saline (c.f., alkaline) lakes found oc- 
casionally throughout this area of North 
Dakota particularly during periods of low 
water. It should be noted that .Alkali Uike 
was not omitted from our enumeration, only 
mis-recognized in terms of size. Our prime 
objective in applying the pro|x>rtion esti- 
mation algoritlim was to recognize many of 
the smaller water features which otherwise 
would not have been detected. From an 
examination of the results, it appears that all 
ponds larger than 0.5 hectare (1.3 acres) 
were consistently recognized and that many 
smaller ponds to 0.13 hectare (0..33 acre) as a 
minimum were recorded. This was a 
threefold improvement relative to the 
single-waveband thresholding technicpie. 
For the 287 km* area illustrated lioth in F'igures 
10 and 1 1, an increase of 131 per cent in the 
numbe* of water bodies tabulated was ob- 
served relative to the niimlrer tabulated by- 
using the thresholding algoritlim. 

Summary and Conclusions 

The mapping of o'len w-ater as an indicator 
of waterfowl habitat quality has been carried 
out by using two different recognition 
techniques, a single waveband thresholding 
approach and a multiple waveband approach 
termed "proportion estimation.” The single 
waveband technique has proven simple to 
implement. Its computer algorithm was 
rapid and accurately recognized prairie lakes 
and large ponds. The resultant products of 


this processing technique were thematic 
maps and statistical tabulations de.scribing 
open surface water conditnns. The maps 
served to portray visually the location and 
fre(|Uency of surface water ho<lies but usu- 
ally necessitated additional interpretation. 
Statistical tabulations provided a useful col- 
lation and data summary. 

The proportion estimation technique, 
which utilizes the increased information 
content of multiple spectral hands, estimates 
the fr.iction of a resolution cell which may he 
composed of open water. This technique al- 
losved for the recognition of a greater 
number of small pimds not previously iden- 
tifu'd iuid improsed tlie apparent spatial res- 
olution of the data by a factor of three in 
comparison to the single w> .eband 
thresholding technicpie. 

Because of the effectiveness of the several 
recognition techniques described, develop- 
ment of an operational system using high- 
altitude sensors for synoptic-ally characteriz- 
ing w-etland habitat conditions appears to be 
a realistic goal. In future operational systems 
which recpiire the delineation of water fea- 
tures smaller than a resolution cell in size, 
we feel that a combination of the single 
w-aveband and the proportion estimation 
techniipies is likely to Iw employed for the 
sake of accuracy and efficiency. Such a dual 
processing approach w-ould elTectively use 
the single w.ivehand technique to identify 
those resolution cells comprised totally of 
w-ater. The proportion estimation algorithm 
which is slow-er running and consequently 
more costly w-ould be used to delineate only- 
water peripheral to already- identified |)onds 
and lakes and those w-ater features w hich are 
less than a resolution cell in size. 
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